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Abstract 



The knowledge of the nuclear symmetry energy of hot neutron-rich matter is important for understanding 
the dynamical evolution of massive stars and the supernova explosion mechanisms. In particular, the electron 
capture rate on nuclei and/or free protons in presupernova explosions is especially sensitive to the symmetry 
energy at finite temperature. In view of the above, in the present work we calculate the symmetry energy as 
a function of the temperature for various values of the baryon density, by applying a momentum-dependent 
effective interaction. In addition to a previous work, the thermal effects are studied separately both in the 
' kinetic part and the interaction part of the symmetry energy. We focus also on the calculations of the mean field 

(-H ' potential, employed extensively in heavy ion reaction research, both for nuclear and pure neutron matter. The 

proton fraction and the electron chemical potential, which are crucial quantities for representing the thermal 
evolution of supernova and neutron stars, are calculated for various values of the temperature. Finally, we 
I construct a temperature dependent equation of state of /3-stable nuclear matter, the basic ingredient for the 

. evaluation of the neutron star properties. 

^ 1 Introduction 



The determination of the nuclear symmetry energy (SE) based on microscopic and/or phenomenological approaches 
is of great interest in nuclear physics as well as in nuclear astrophysics. For instance, it is important for the study of 
, the structure and reactions of neutron-rich nuclei, the Type II supernova explosions, neutron-star mergers and the 
stability of neutron stars. In addition, the SE is the basic ingredient for the determination of the proton fraction 
and electron chemical potential. The above quantities determine the cooling rate and neutrino emission flux of 
■ protoneutron stars and the possibility of kaon condensation in dense matter [U . 

Heavy-ion reactions are a unique means to produce in terrestrial laboratories hot neutron-rich matter similar 
to those existing in many astrophysical situations [3]. Although the behavior of the SE for densities below the 
saturation point still remains unknown, significant progress has been made only most recently in constraining the 
SE at subnormal densities and around the normal density from the isospin diffusion data in heavy-ion collisions 
[H [S] . This has led to a significantly more refined constraint on neutron-skin thickness of heavy nuclei [51 [7] and 
the mass-radius correlation of neutron stars fS] . For densities above the saturation point the trend of the SE is 
model dependent and exhibits completely different behavior. 

Up to now the main part of the calculations concerning the density dependence of the SE is related with the cold 
nuclear matter (T = 0). However, recently, there is an increasing interest for the study of the SE and the properties 
of neutron stars at finite temperature [21 M O El dH [HI [SI [E] ■ The motivation of the present work is to clarify 
the effects of finite temperature on SE and to find also the appropriate relations describing that effect. Especially 
we focus on the interaction part of the SE, where so far it has received little theoretical attention concerning its 
dependence on the temperature. 

In order to investigate the thermal properties of the SE, we apply a momentum dependent effective interaction 
model. In that way, we are able to study simultaneously thermal effects not only on the kinetic part of the symmetry 
energy but also on the interaction part. The present model has been introduced by Gale et al. [T^l [13 [TSl [TH] in 
order to examine the influence of momentum-dependent interactions on the momentum flow of heavy ion collisions. 
Over the years the model has been extensively applied in the study not only of the heavy ion collisions but also in 
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the properties of nuclear matter by a proper modification [20l [211 [111 [23] . A review analysis of the present model 
is presented in Refs. [21 [18] . 

In the present work we study the thermal properties of the nuclear symmetry energy by applying the above 
phenomenological model focusing mainly on the temperature dependence of the kinetic and interaction part of the 
SE as well as the total SE. Though it is well known how the temperature affects the kinetic part of the symmetry 
energy [31 [^H [53] the temperature dependence of the interaction part of the SE has so far received little theoretical 
attention. In addition, we determine the temperature dependence of the proton fraction as well as of the electron 
chemical potential. Both of the above quantities are related with the thermal evaluation of the supernova and the 
proton-neutron stars. The single particle potential for the pure neutron matter and the symmetric nuclear matter, 
extensively applied in heavy ion collision research, is also estimated for various values of the temperature. Finally, 
we construct the equation of state (EOS) of /3-stable matter which is the basic ingredient for calculations of the 
neutron star properties. 

The plan of the paper is as follows. In Sec. II the model and the relative formulae are discussed and analyzed. 
Results are reported and discussed in Sec. Ill, while the summary of the work is given in Sec. IV. 



2 The model 

The schematic potential model, used in the present work, is designed to reproduce the results of the more mi- 
croscopic calculations of both nuclear and neutron-rich matter at zero temperature and can be extended to finite 
temperature [2j. The energy density of the asymmetric nuclear matter (ANM) is given by the relation 



e{nn,np,T) = ^ktni^n, T) + el^„{np,T) + Vint{nn,np,T), 



(1) 



where n„ (rip) is the neutron (proton) density and the total baryon density is 
the kinetic parts are 



n — Un + Tip. 



e'^,„{n,„T) + elJnp,T) = 2 



(27r)3 2m 



(/«(nn, k,T) + fp{np,k,T)), 



The contribution of 



(2) 



where /,-, (for t = n,p) is the Fermi-Dirac distribution function with the form 



fr{nr,k,T) = 



1 -I- exp 



3i-(n, fc,r) - jjLr(n,T) 



(3) 



The nucleon density tIt is evaluated from the following integral 



1 + exp 



er{n,k,T) - iir{n,T) 



(4) 



In Eq. ([3]), er{n,,k,T) is the single particle energy (SPE) and ii.r{n,T) stands for the chemical potential of each 
species. The SPE has the form 

erin,k,T) = ^— + Ur{n,k,T), (5) 
zm 

where the single particle potential Ur{'n,k^T), is obtained by differentiating Vint i-e. Ut = dVintinn,iip,T)/dnT- 
Including the effect of finite-range forces between nucleons, in order to avoid acausal behavior at high densities, 
the potential contribution is parameterized as follows [2] 



Vint{nn,np,T) = -^Ano 
2 ^ 



iBno [|-(l + a;3)(l-2x)^] u^+^ 
1 + lB' [|-(i+a;3)(l-2a;)2] 



d^k 



{2C,+AZi) 2 / ^g(fc,A,)(/„ + /p) 



/Ail. 
^^g{k,Ki){Ul-x)+fpX) 



(6) 
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where x = Up/n is the proton fraction and u = n/riQ, with no denoting the equihbrium symmetric nuclear matter 
density no = 0.16 fm~^. The constants A, B, a, Ci, C2 and B', which enter in the description of symmetric nuclear 
matter and the additional parameters xo, X3, Zi, and Z2 , used to determine the properties of asymmetric nuclear 
matter, are treated as parameters constrained by empirical knowledge [2]. The function g(k,Ai) suitably chosen 
to simulate finite range effects is of the following form 



(7) 



where the finite range parameters are Ai = l.S/c^ and A2 ~ 3kp and k^ is the Fermi momentum at the saturation 
point no- 

The entropy density Sr{n, T) required for the calculations of the total pressure and for the EOS, has the same 
functional form as that of a non interacting gas system, that is 

/d^k 
[fr In U + (1 - fr) ln(l - /.)] . (8) 

The ratio entropy/baryon is given by Sr(n,T) — ST{n,T)/n. The baryon pressure Pb{n,T), needed to construct 
the EOS, is given by 

Pb{n,T)^T ^ Sr{n,T) + ^ n^/x^(n, T) - ea„„i(n, T). (9) 

T—p.n T—p,n 

Finally, the total energy density and pressure of charge neutral and chemically equilibrium nuclear matter are 

etotin,T) = ei,in,T) + ^ eiin,T), (10) 

l—e~ 

Ptot{n,T) = Pbin.T) + ^ Pi{n,T). (11) 

l—e~ ,fi~ 

The leptons (electrons and muons) originating from the condition of the /3-stable matter are considered as non- 
interacting Fermi gases. 

The above analysis holds in general for the asymmetric nuclear matter. Below, in order to calculate the thermal 
effect on the SE, we will focus our study on two cases, i.e. the symmetric nuclear matter (SNM) and the pure 
neutron matter (PNM). 

2.1 Symmetric nuclear matter 

The energy density of SNM is given by Eqs. ([T]) and ^ by setting x = 1/2, that is [2] 

f (fk h^k^ f (fik ri^k^ 1 ^ 2 Bnou'^+i 

es„m(n,T) = 2 / -pr-r^^ — Jn + 2 / jtt^^ — fp + t^^t^ou ' 



u 



J (27r)3 2m ' J (27r)3 2to -"^ 2 ^ l + B'u"-^ 

J^gik, A0/« + uJ2a2 j—^gik, AO/p. (12) 



In addition, the single particle potential C/J„,„(n, fc, T) in the case of SNM, defined from the relation C/J,„„ = 
dVsnm/dnr, is easily calculated and given by 



i=l,2 



A'' 

A,: 



(13) 



It is obvious from Eq. ([T^ that C/J„„j(n, k, T) is separated in two terms. The first one corresponds to the momentum 
independent part, while the second one corresponds to the momentum dependent one. The term UJ„j^{n,T) has 
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the following form 



UIn,nin,T) ^ Au + 



Bu^ia + l + 2B'u''-^) 



(1 +B'U'^-1)2 



no J 27r 



fr. 



T — p,n. 



(14) 



At zero temperature {T — 0), where 6{kp^ — k), the integrals in Eqs. and are calculated analytically 
(see Appendix A for more details). 

2.2 Pure nuclear matter 

The energy density of PNM is given by Eqs. ^ and ([6]) by setting a; = and fp — 0, that is [2] 



d^k h\\ 1 ^ ,2 |Sno(l-X3)M"+i 

(27r)'^ 2m 3 1 + ^i?'(l - X3)u'^-^ 



(2^ 



g{k,A,)f„. 



(15) 



The single particle potential Upnj^{n,k,T) in the case of PNM is defined from the relation J7p„„ — dVpum/dun is 
written as 



C/;„„(n, fc, T) = [7;_(n, T) + ^ (3G - 4Z,) 



(16) 



The momentum-independent part is 



-A{l-xo)u + 



B{\ - a;3)M° 



[l + fB'(l-X3)w'^-l]2 



{a + l) + -B'{l-x^)u--^ 



(17) 



The integrals in Eqs. (fT51) and ([T7]), similarly to the case of SNM, at T = are calculated analytically (see 
Appendix A for more details). 

2.3 Asymmetric nuclear matter-Nuclear symmetry energy 

The energy density of ANM at density n and temperature T, in a good approximation, is expressed as 

eanm{n,T,x) = es„„(n,T,x = 1/2) + esy„iin,T,x), (18) 

where 



in,T,x) = nil - 2xfElllin,T) = n{l - 2xf {E%l,{n,T) + ETyl{n,T)) 



(19) 



In Eq. (fT9)) the nuclear symmetry energy El°l^{n,T) is separated in two parts corresponding to the kinetic con- 
tribution E^y"^ (n, T) and the interaction contribution E"™^ (n, T) . In the present work we will concentrate on the 
systematic study of the thermal properties of the above two quantities. 

From Eqs. p8|) and (fT9)) and setting a; = we obtain that the nuclear symmetry energy i?*°^(n, T) is given by 



Elylnin, T) = - {Cpnmin, T) - es„m(r7,, T)) . 



(20) 



Thus, from Eqs. (jl2p and (jlSp and by a suitable choice of the parameters xq, X3, Z\ and Z2, we can obtain different 
forms for the density dependence of the symmetry energy E^°y^{n,T^. It is well known that the need to explore 
different forms for E^°y^(n^ T) stems from the uncertain behavior at high density [2]. In the present work, since we 
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are interested mainly in the study of thermal effects on the SE, we choose a specific form of the SE enabling us to 
reproduce accurately the results of many other theoretical studies [26]. According to this choice the SE, at T = 0, 
is expressed as 




Kinetic Interaction Kinetic Interaction 



where the contributions of the kinetic and the interaction term are separated clearly. The parameters xq, X3, 
Z\ and Zi are chosen in order that Eq. (PO)) . for T = 0, to reproduce the results of Eq. (PT|) . In addition, the 
parameters B, a, Ci, C2 and B' are determined in order that E{n = no) — mc^ = —16 MeV, uq = 0.16 fm^^, 
and the incompressibility to be Kq = 240 MeV. 

The single particle potential U^^^{n,k,T), in the case of ANM defined from the relation U^^^^ — dVanm/driT, 
is written as 

dV 

U:„^{n, k, T) = C/;„„(n, fc, T) + = UJ^,^in, k, T) + UJy^{n, T, x), (22) 

where 

V,ym{n,T,x) - (1 - 2x)'nETyl{n,T). (23) 

It is easy to find that the term Uly^{n,T), in the case of T = and by applying expression (PT|) . is given by (see 
alsoref. [27]) 

C/J^„(n,T,x)=±34u(l-2x), (24) 

where + and— stand for neutrons and protons respectively. In the general case where thermal effects are included 
in our calculations, the i?™^(n, T) takes the form 

E^:il{n,T) ^ au\ (25) 

where a and 6 are temperature dependent constants (see Eq. (|¥T]l on Sec. III). Thus, after some algebra, we get in 
a good approximation, the relation 

Uly^{n,T,x) ~ ±2au''(l - 2x). (26) 

The above relation is needed for the calculation of the single particle energy er{n^ k, T) in the /3-stable matter and 
afterwards for the calculation of the Fermi-Dirac function fr{n,T) which is the basic ingredient for the determina- 
tion of the entropy density Srin, T). 



2.4 Proton fraction-Electron chemical potential 

The key quantity for the determination of the equation of state in /3-stable matter is the proton fraction x, which 
is a basic ingredient of Eq. . In /3-stable matter the processes [55] 

n — >p + e^ + De, p + e~ — >n + iye, (27) 

take place simultaneously. We assume that neutrinos generated in these reactions have left the system. This implies 
that 

fj. ^ fj.n - ^J.p = fJ-e, (28) 

where ^m^^p and jjte are the chemical potentials of the neutron, proton and electron respectively. Given the total 
energy density e = e(n„,np), the neutron and proton chemical potentials can be defined as 

de de 

Hence we can show that 

de/n BE 

H = lj,n-fj,p = — _|„ = __|^. (30) 

In (3 equilibrium one has 

f)E B 

-g^ = -g^{Eb{n,x)+E,{x))^0, (31) 
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where Eh{n^ x) the energy per baryon and Ee{x) the electron energy. The charge condition impUes that rig 
nx 01 kp^ = kpp- Combining relations (fT8|) . (fT9|) and ([30|) we get 

Me(n, T) = fi{n, T) = 4(1 - 2x)Etl,{n, T). 



(32) 



From Eq. (j32p it is obvious that the proton fraction a; is not only a function of the baryon density n but, in addition, 
depends on the temperature T i.e. x = x{n,T). 

For relativistic non-degenerate free electrons we have 



Tie = xn 



(27r3) 



1 + exp 



T 



(33) 



Or, using Eq. (|32p and performing the angular integration we get 



1 



k^dk 



1 + exp 



V^^fe^c^+mfc* - 4 ( 1 - 2 a; ) _E|2 m ( " ' "2") 



(34) 



Eq. p4p determines the equilibrium electron (proton) fraction x{n,T) since the density and momentum dependent 
symmetry energy El°l^{n,T) is known. 



2.5 Calculations recipe 

We focus our attention on the calculation of the El°l^{n,T) with the help of Eq. (120]) . Thus, one has to calculate 
first the energy densities in pure and in symmetric nuclear matter as a function of the density n and for fixed values 
of temperature T. As an example of the calculations procedure at finite temperature (the results for T = are 
included in the Appendix A), we consider the case of pure neutron matter. The procedure is similar in the case of 
symmetric nuclear matter (see Ref. [2 ). 

The outline of our approach is the following: For a fixed neutron density n„ and temperature T, Eq. (|4]) may 
be solved iteratively in order to calculate the variable 



ri{n;T) 



firjn-yT) - U{n\T) 
T 



(35) 



The knowledge of 'q{n, T) allows the last term in Eq. p7|) to be evaluated, yielding f7(n; T) which may then be 
used to infer the chemical potential from 



lJir{n;T) = Tf^{n-T) + U{7i-T), 



(36) 



required as an input to the calculation of the single particle spectrum erin, k, T) in Eq. ([5]). Using erin, k] T), the 
energy density in Eq. (fT5)) is evaluated. 



3 Results and Discussion 

According to our calculation recipe, given in the previous subsection, we calculate the energy densities of PNM 
and SNM as functions of the density, for various values of the temperature T. As a second step, we calculate the 
Elyl^{n,T) from Eq. d^Ol)- The knowledge of El°*^{n,T) is required for the evaluation of the proton fraction x 
from Eq. as well as for the electron chemical potential fie = fi from Eq. ([5^ . Finally from Eqs. (O, PTI)) and 
([TT|) we construct the EOS of /3-stable matter for various values of the temperature T. It is worth pointing out 
that in the present work we do not include the muon case, since we restrict ourselves mainly on the temperature 
dependent behavior of the SE. According to our plan, in future work we will extend the treatment to include also 
the muon case in order to study the detailed composition and the thermal properties of neutron-rich matter with 
applications in neutron star structure and thermal evaluation. 

In Fig. 1 we check the validity of approximation p^ . We plot the difference E{n,T,x) — E{n,T,x — 1/2) as 
a function of (1 — 2x)^ at temperature T = 0, T = 20 and T = 50 MeV for three baryon number fractions i.e. 
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u = 1, u = 2 and u = 3. It is seen that an almost linear relation holds between E{n, T, x) — E{n, T,x = 1/2) and 
(1 — 2a;)^, even closer to the case of pure neutron matter {x = 0), indicating the validity of approximation (fT8|) . 

In Fig. 2 we indicate the behavior of the SE as a function of the temperature T for various fixed values of the 
baryon density n. More precisely, in any case, we plot i?*°^(T;n), as well as £'^*^(T;n) and -E™^(T;n) as a 
function of T for n = 0.1,0.2,0.3,0.5 fm~^. The most striking feature of the above analysis is a decrease of the 
SE (total, kinetic and interaction part) by increasing the temperature. This is consistent with the predictions of 
microscopic and/or phenomenological theories [3l IT3l [14] 

In order to illustrate further the dependence of the symmetry energy on the temperature and to find the 
quantitative characteristic on this dependence, the values of Esym{T\n) for various values of the density n are 
derived with the least-squares fit method and found to take the general form 



A 

1 + (TVTo)- 



Esy.n{T; n) = ^ , ^, + B. (37) 



The values of the density dependent parameters A, B, Tq and c, for El°y^{T;n), E^y^{T]n) and Elyl^{T;n) 
for n = 0.1,0.3,0.5 fm~^ are presented in Table 1. It is easy to find that in the case of low temperature limit 
{T/Tq <C 1) all kinds of the symmetry energy decrease approximately according to Esym (T; n) cx Ci — C2T^ (where 
Ci and C2 density dependent constants). In the high density limit {T/Tq ^ 1) the symmetry energy decreases 
approximately according to Esym{T;n) oc CaT"^ + C4 (where also C3 and C4 are density dependents constants). 
It is noted that the same behavior holds for El°l^(T] n) as well as for E^y^{T; n) and E™l^{T; n). This behavior is 
well expected for the kinetic part of the symmetry energy (see also Ref. [SJIl^), where analytical calculations are 
possible (see the prove in Appendix B). From the above study, it is concluded that there is a similar temperature 
dependence both for the kinetic and the interaction part of the symmetry energy and consequently for the total 
symmetry energy, in the case of momentum dependent interaction. Recently, the temperature dependence of the 
kinetic and interaction part of the SE has been studied and illustrated in Ref. [13] ■ The results of the present work 
agree with those of Ref. [Tl] although different models have been employed to evaluate SE. 

In Fig. 3, we plot £^*°^(T;n) as a function of temperature for various low values of the baryon density. In 
the same figure we also include experimental data of the measured temperature dependent symmetry energy from 
Texas A&M University (TAMU)[29] and the INDRA-ALADIN Collaboration at GSI [30|. The comparison then 
allows to estimate the required density of the fragment-emitting of the experiments. As pointed out by Li et.al. [3] 
the experimentally observed evolution of the SE is mainly due to the change in density rather than temperature. 

Fig. 4 illustrates the behavior of the i^*°'^'(n; T) (a), E^}^^^{n; T) (b), Ell^*^{n; T) (c), as a function of the baryon 
density n for various fixed values of the temperature T. The case T = corresponds to the fundamental expression 
of the present work i.e. 

i;*°*„(u; T = 0) = 13^2/3 + i7y. (38) 

In any case, the trends of the various parts of the symmetry energy are similar. An increase in the temperature 
leads just to a shift to lower values for the symmetry energy. It is worth pointing out that, the maximum decrease 
oiEll\^{n; T), in the area under study (for T = McV up to T = 50 McV), is between 40% (for n = 0.1 fm'^) and 
4% (for n — Ifm^'^). Correspondingly, the decrease of E^y^{n;T) is between 57% (for n = 0.1 fm~^) and 5% (for 
n = l fm-3) and of the E'J^l^in; T) is between 22% (for n = 0.1 fm^^) and 5% (for n = 1 fm'^). It is obvious that 
the thermal effects are more pronounced on the kinetic part than in the interaction part of the symmetry energy 
and in addition, more pronounced in lower values of the baryon density. 

The total symmetry energy Ely\^(u;T), for various values of the temperature T, was derived with the least- 
squares fit on the numerical results taken from Eq. (|20p and has the form 



^(u;r = 5) = 1.676 + 29.711m- 2. 110?i2 + 0.275w^-0.015u'*, 
El°yl^{u;T ^ 10) = -0.118 + 30.863m - 2.455^2 + 0.3257/3 -0.017^^^ 
Esym{u;T ^ 20) = -1.910 + 29.470m - 1.466m^ + 0.120m3-0.004u'^ 



sym \ 

ot . 

sym \ 

sym ( 

(m;T = 50) = 0.099+ 18.172m +2.9m2- 0.548m3 + 0.033m^ (39) 



sym \ 

It is also useful to record some relations for Elyl^{u; T) derived by least-squares fit on the numerical results, in 
the case where SE is parameterized in a way similar to that one holding for T = 0. In that case, the parametrization 
is the following (the case E\'^^{u\ T = 0) is included also for comparison) 

E%Uu;T = 0) = 13m2/3 + 17m, 
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rptot (^,. 'T' ~ 
^syraK^-! ~ 


= 5) 


rptot 7^ 

^sym\^^ ^ ~ 


0) 


- 0.374 


u' 


-0.956 


rptot . 7^ 


10) 


rptot (^.. 'T' 

^symK^T-^ ~ 


0) 


- 1.235 


u 


-0.804 




20) 


rptot (^,. 'T' 

~ ^symK^-! ^ ~^ 


0) 


- 3.420 


u 


-0.520 


-rptot . 7^ 


50) 


rptot (^,. 'T' 

^symK^-! ^ ~ 


0) 


- 9.300 


u 


-0.097 



(40) 

From Eq. ([301), the decrease of the SE as a result of increasing T, is evident. 

The interaction part of the symmetry energy i?™*„(u;r) for various values of the temperature T was derived 
by a least-squares fit on the numerical results taken from Eqs. and (PO)) and has the form 

ETyL{u-T = b) = 17.041 
i;™^(u;T=10) = 16.782 ui "05^ 
£»^(m;T = 20) = 16.022 

£;™^(u;T = 50) = 13.404 M^ iO'*. (41) 

Similarly, for the kinetic part of the symmetry energy £'g™j(it;T) we obtain 

E%l^{u;T = b) = 12.856 ^t°•^^^ 

£;,*^^^(u;T= 10) = 12.504 

E%^^{u;T ^ 20) = 11.518 7.°-^36^ 

E^;i,{u;T = 50) = 8.577 (42) 

In Fig. 5 we plot the total energy per particle of the PNM (a) and of the SNM as a fmiction of the density for 
various values of the temperature. In both cases it is concluded that the thermal effects become more pronounced 
when T > 10 MeV and for baryon densities n < 0.5 fm~'^. 

Fig. 6 displays the single particle potential Upnm{n,T, k) of the PNM as a function of the momentum k for 
various values of the density n and temperature T. An increase of T leads to corresponding increase of the values of 
the Upnm{n, T, fc), an effect, expected to be more pronounced for lower values of the baryon density (n = 0.1 fm~'^) 
compared to highest (n — 0.5 fm^^). The same trend holds also for the single particle potential Usnm{n,T, k) of 
the SNM plotted in Fig. 7. Observing Figs. 6 and 7 one might expect that the change of T will affect slightly the 
nucleons with high momentum k. This could be seen by plotting the single particle energy e^-(f^, fc, T) (see Eq. ([5])) 
as a function of k. However, the above effect cannot be seen in the present work, where we plot just the single 
particle potential U'^{n, fc, T) as a function of k. 

In Fig. 8 we display the single particle potential of neutron U'"'{n,T, k) (Fig. (a),(b)) and proton UP{n,T, k) 
(Fig. (c),(d)), in /3-stable matter, as a function of the momentum k for various values of the temperature T for 
n = 0.1 and n = 0.5 fm~^. The potential U"^{n,T, k) is evaluated according to Eq. ([22l The most striking feature 
of Fig. 8 is the reduced thermal effect for high values of the baryon density, especially in the case of the neutron 
single particle potential. In the case of the proton, thermal effects are more pronounced. 

In Fig. 9(a) the proton fraction x is displayed, calculated from Eq. p4)) as a function of n for various values of 
T. Thermal effects increase the value of x between 57% (for n = 0.1 fm~'^) and 2% (for n — 1 fm~^). This effect 
is directly related with the dependence of x on the symmetry energy. As discussed previously, the temperature 
influences slightly the symmetry energy at high values of the density and consequently this is reflected in the values 
of X. It is stressed that x depends on T in two ways, as one can see from Eq. ((M)) . That is, it depends directly on T 
due the Dirac- Fermi distribution and also depends on the symmetry energy which is also temperature dependent. 

In Fig. 9(b) we present the electron chemical potential /ie as a function of the density n for various T. An 
increase of T decreases /ig. The effect is more pronounced when T > 20 MeV. We mention that the rate of electron 
capture on both free and bound protons depends in a very sensitive way on the difference /} = /i„ — /ip = /ig 
between neutron and proton chemical potentials [5]. Larger values oi jl = fie inhibit the neutronization process, 
since it becomes more difficult to transform a proton into a neutron. 

Finally, in Fig. 10 we present the equation of state of beta stable matter constructed by applying the present 
momentum-dependent interaction model, for various values of the temperature T. It is obvious that the thermal 
effects are enhanced when T > 20 MeV. The above EOS is very important for the calculation of the neutron stars 
properties and also in combination with the calculated proton fraction and electron chemical potentials for the 
thermal evaluation of the neutron stars. 
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4 Summary 



The knowledge of the nuclear symmetry energy of hot neutron-rich matter is important for understanding the 
dynamical evolution of massive stars and the supernova explosion mechanisms. In view of the above statement, we 
investigate, in the present work, the thermal effects on the nuclear symmetry energy. In order to perform the above 
investigation we apply a model with a momentum-dependent effective interaction. In that way, we are able to study 
the thermal effect not only on the kinetic part of the symmetry energy but also on the interaction part which, in 
turn, due to a momentum dependence, is affected by the variation of the temperature. It is concluded that, in 
general, by increasing T we obtain a decreasing SE. Our finding that both kinetic and interaction parts exhibit the 
same trend both for low and high values of the temperature is an interesting result. Analytical relations, derived 
by the method of least squares fit are given also for the above quantities. Temperature effects on the pure neutron 
matter and also on symmetric nuclear matter are also investigated and presented. The single particle potential of 
proton and neutron is of interest in heavy ions collisions experiments, is calculated also for pure neutron matter, 
symmetric nuclear matter and /3-stable matter for various values of the baryon density and fixed values of T. It 
is concluded that thermal effects are more pronounced for low values of the density n, where for high values of n 
the effects are almost negligible. Quantities, which are of great interest for the thermal evaluation of supernova 
and neutron stars, i.e. the proton fraction x = x{n,T) and the electron chemical potential /ie — fiein,T), are 
calculated and their temperature and density dependence is investigated. Thermal effects are larger for low values 
of the density and high values of T. 



Appendix A 

The energy density of the SNM as well as of the PNM, at zero temperature are easily calculated from Eqs. (fT2|) 
and Uni) respectively by setting = 0{kp^ — k) (where 6{kF^ — k) is the theta function and kp^ is the Fermi 
momentum of the nucleon r) and takes the following forms 



,(n,fc;r = 0) 



1,2 \ F / '\ h. „ k. r 



5 3 1 + |i3 (1 — X3)u'^~^ 

+ 5„„„i:,3a-4..,(^)7e#!-u„-.e|:^!), (44, 

where E% = U^ky /2 m is the Fermi energy of nuclear matter at the equilibrium density. 



Appendix B 

In order to compare the numerical results obtained from the kinetic part of the symmetry energy E^^^{n^T) with 
those predicted from analytical calculations, we calculate El°l^' {n, T) in the low and in the hight temperature limit 
as follows 

Low temperature limit 

The kinetic energy per nucleon EJ.^^{n,T) at low temperature (T ^ Ep) has the form [311 [321133] 



:E^P 



5 n 1 T 

El 



12 



(45) 
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where Ep = {hkp)'^/2m = h^iSn'^rir)'^^^ /2m. Considering that S = I — 2x = (n„ — np)/{nn + rip) after some 
algebra we found that the Ekin{n, T, 5) of a two-component Fermi gas has the form 

^ ' (|t)'((1 + <5)V3 + (i_ 5)1/3), (46) 



10 {E_ 



F) 



where {Ep) =^ 3/5Ep. Expanding expression (|46p aromid the symmetric point (5 = or a; = 1/2 the kinetic energy 
takes the approximated form 

r) . iE,) + A^r^ + (1 - 2.r- [l(E,) - 1^T=) , (47) 

" v ' 

with the contribution of the symmetry energy written exphcitly. It is obvious that in the low temperature limit 
E%l^{n,T) behaves as E^;"^in,T) ex Ci - CsT^. 



High temperature limit 

The kinetic energy per nucleon EkininjT, S) of a two-component Fermi gas at high temperature (T ^ Ep) is 



replaced by a virial expansion in nX^ where A — y2Trh^/mT is the quantum wavelength. So, £'fci„(n,T) is given 
by the relation |32l [25] 



E,.Un,S,T) = -T^ a ( — j ((1 - S)'^+' + (1 + S)'^+') . (48) 

Expanding expression (|48p around the symmetric point (5 = Oora; = l/2 the kinetic energy takes the approximated 
form 



EkUn,T,S) = 



It is seen that in the high temperature hmit E^l'^{n,T) behaves as E^l'^{n,T) cx CiT'^/'^ + C2T-'^ H . 
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Figure 1: The difference E{n, T, x) — E{n, T,x — 1/2) as a function of (1 — 2a;)^ at temperature T = 0, T = 20 and 
T = 50 MeV for three baryon number fractions u — 1, u — 2 and u — 3. 



Table 1: The values of the density dependent parameters A, B, Tq and c, for El°l^{u; T), E';il^{u- T) and Eif^{u] T) 
for n = 0.1, 0.3, 0.5 fm~^. For more details see text. 







n^O.l fm-^ 






n=0.3 fm"^ 






n=0.5 fm 




Parameters 


rptot 
sym 


rpkin 
sym 


rpint 
sym 


rptot 
sym 


rpkin 
sym 


rpi7it 
sym 


rptot 
sym 


sym 


rpint 
sym 


A 


10.105 


7.864 


2.559 


17.079 


15.230 


10.832 


44.164 


20.548 


19.442 


B 


9.969 


1.679 


7.969 


19.328 


4.504 


21.240 


36.895 


7.162 


33.887 


To 


25.692 


30.549 


19.027 


41.004 


73.143 


47.772 


57.011 


109.193 


73.551 


c 


1.610 


1.518 


1.866 


1.856 


1.904 


1.992 


1.982 


2.156 


2.026 
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Figure 2: Temperature dependence of the total nuclear symmetry energy and its interaction and kinetic energy 
part for various values of the baryon density n. 
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Figure 3: Temperature dependence of the symmetry energy for low values of the baryon density (n = 
0.06,0.08,0.10,0.12,0.14 fm-^). The experimental data are from Ref. [55] (solid squares) and Ref. (open 
squares) are included for comparison. 
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n(fm'') n(fni') n(fm"') 



Figure 4: Density dependence of the total nuclear symmetry energy i?*°^(n,T) as well its kinetic E^^{n,T) and 
interaction E™l^{n,T) part for various values of the temperature T. 
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Figure 5: (a) The energy per particle of pure neutron matter as a function of the baryon density for various vakies 
of the temperature T. (b) The energy per particle of symmetric nuclear matter as a function of the baryon density 
for various values of the temperature T. 
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Figure 9: (a) The proton fraction x in /^-stable matter as a function of the density n for various values of the 
temperature T. The straight Une corresponds to the case x = 11%. (b) The electron chemical potential /ie = A = 
— /Ltp as a function of the density n for various values of the temperature T. 




Figure 10: The equation of state P = P{e) of /3-stable matter corresponding to the present momentum-dependent 
effective interaction model for various values of the temperature T. 
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